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Abstract. We report on an approximately twelve hour long X-ray flare from the low-mass X-ray binary 
KS 1731—260. The flare has a rise time of less than 13min and declines exponentially with a decay time of 
2.7 hours. The flare emission is well described by black-body radiation with peak temperature of 2.4 keV. The 
total energy release from the event is 10 42 erg (for an assumed distance of 7kpc). The flare has all the charac- 
teristics of thermo- nuclear X-ray bursts (so-called type I X-ray bursts), except for its very long duration and 
therefore large energy release (factor of 1500-4000 longer and 250-425 more energy than normal type I X-ray 
bursts from this source). The flare is preceded by a short and weak X-ray burst, possibly of type I. Days to weeks 
before the flare, type I X-ray bursts were seen at a rate of ~3 per day. However, after the flare type I X-ray 
bursting ceased for at least a month, suggesting that the superburst affected the type I bursting behaviour. The 
persistent emission is not significantly different during the non-bursting period. We compare the characteristics of 
this event with similar long X-ray flares, so-called superbursts, seen in other sources (4U 1735—44, 4U 1820—30, 
4U 1636-53, Ser X-l, GX 3+1). The event seen from KS 1731-260 is the longest reported so far. We discuss 
two possible mechanisms that might cause these superbursts, unstable carbon burning (as proposed recently) and 
electron capture by protons with subsequent capture of the resulting neutrons by heavy nuclei. 

Key words, accretion, accretion disks — binaries: close — stars: individual (KS 1731— 260) — stars: neutron — 
X-rays: bursts 



1. Introduction 

Many low-mass X-ray binaries show thermo-nuclear ex- 
plosions, or type I X-ray bursts (for a review see Lewin et 
al. 1993). These appear as rapid (^1 s) increases in the X- 
ray flux, followed by an exponential decline, with typical 
durations of the order of seconds to minutes. The (net) 
burst spectra are well described by black-body emission 
from a compact object with ~10km radius and tempera- 
ture of ~l-2keV. The inferred temperature decreases dur- 
ing the decay, indicating cooling of the neutron star sur- 
face. Typical integrated burst energies are in the ~10 39 to 
10 40 erg range. 

Recently, new modes of unstable burning on neutron 
stars may have been revealed, with the discovery of X- 
ray flares lasting several hours (Cornelisse et al. 2000; 
Strohmayer 2000; Strohmayer & Brown 2001; Wijnands 
2001; Kuulkers 2001; Cornelisse et al. 2001). The re- 
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currence times of these events are not well constrained. 
Wijnands (2001) reported two events from 4U 1636—53 
which were ~4.7 years apart. So far, only seven events have 
been found, indicating that they are rare. The availability 
of X-ray instruments such as the BeppoSAX/WFC and 
the RXTE/ASM, which more frequently monitor the X- 
ray sky, and the RXTE/PCA, which performs long studies 
of the X-ray burster population, is the reason that these 
events have now been discovered. 



These long X-ray flares share many of the charac- 
teristics of type I bursts, which is the reason they are 
attributed to thermonuclear runaway events (see e.g. 
Cornelisse et al. 2000). What distinguishes them from nor- 
mal type I bursts are the long duration (exponential de- 
cay times of a few hours), the large fluences (~10 42 erg), 
and the aforementioned extreme rarity. Because of the 
large fluences, they are also referred to as 'superbursts' 
So far, these events have only been observed in sources 
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with persistent pre-burst luminosities of ~0. 1-0.3 times 
the Eddington luminosity (e.g. Wijnands 2001). 

The current view is that the superbursts are caused 
by the unstable burning of the ashes of the (un)stable 
hydrogen and/or helium burning (Cumming & Bildsten 
2001; Strohmayer & Brown 2001). Such bursts in principle 
thus not only tell us about properties of material buried 
below the hydrogen and/or helium layer, but also about 
the burning of the hydrogen and/or helium layer itself. 

We here report on a superburst seen with 
the BeppoSAX Wide Field Camera (WFC) from 
KS 1731-260 on 1996 September 23. In Sect. 2, we 
briefly review what is known about KS 1731—260. Next, 
we describe the observations (Sect. 3) and the results 
(Sect. 4). We conclude with a discussion of possible 
mechanisms for the superbursts. Hereafter we use the 
term (X-ray) bursts, whenever we refer to type I (X-ray) 
bursts (unless it is needed for clarity). 

2. KS 1731-260 

KS 1731—260 was discovered as a transient and burst- 
ing source in 1989 by COMIS/TTM onboard Mir-Kvant 
(Sunyaev et al. 1989, 1990). A possible infrared counter- 
part was identified recently using a precise localisation by 
Chandra (Wijnands et al. 2001b; Revnivtsev & Sunyaev 
2002; Groot et al. 2001, in preparation). The source is an 
X-ray burster (Sunyaev et al. 1990; Muno et al. 2000), 
which identifies the compact object as a neutron star. 
Highly coherent X-ray oscillations were discovered by the 
RXTE/PCA at a period of 1.9 ms during several X-ray 
bursts (Smith et al. 1997; Muno et al. 2000) as well as 
kHz quasi-periodic oscillations at ^900 Hz and ~ 1200 Hz 
in the persistent emission (Wijnands & van der Klis 1997). 
The distance to the source, as derived from the peak lumi- 
nosity observed in radius-expansion bursts, is about 7kpc 
(Muno et al. 2000). 

Fig. |l|a shows the RXTE/ASM light curve of 
KS 1731-260 from near the start of the RXTE mission, 
i.e. 1996 Jan 6, to 2001 May 1. The source is highly vari- 
able on time scales of a day to weeks. During the first two 
and a half year of the RXTE mission the overall source 
intensity varied slowly on time scales longer than a week. 
Subsequently, it showed a two and a half year interval 
with strong flares or outbursts on time scales of weeks to 
a month, reflecting a change in source state (see also Muno 
et al. 2000; Revnivtsev & Sunyaev 2001). At the beginning 
of 2001, the source became quiescent (see Wijnands et al. 
2001a). 

3. Observations 

BeppoSAX was launched on 1996 April 30. Onboard are 
two Wide Field Camera's (WFCs); they consist of coded 
mask aperture cameras (Jager et al. 1997). They point 
in opposite directions of one another and perpendicular 
to the Narrow-Field Instruments on the same satellite 
(Boella et al. 1997). The field of view of each WFC is 



40x40 square degrees full- width to zero response with an 
angular resolution of about 5 arc minutes in each direc- 
tion. The band pass of the instruments is 2 to 28keV. The 
imaging capability and sensitivity allow an accurate mon- 
itoring of complex sky regions, like the Galactic Centre. 
Two times per year, in spring and autumn, for about 
two months a dedicated program is carried out to mon- 
itor the Galactic Centre region typically for about half 
a day on a weekly basis. During the performance veri- 
fication phase, however, the Galactic Centre was moni- 
tored continuously (only interrupted by South Atlantic 
Anomaly passages and earth occultations) for an interval 
of ^9 days, four months after launch (see e.g. Cornelisse et 
al. 2000; Kuulkers et al. 2000). KS 1731-260 lies close to 
the Galactic Centre (I = 1.07°, b = +3.66°), ensuring the 
source is well covered. In this paper we use the data from 
from either WFC unit 1 or 2, from the Galactic Centre 
monitoring program as well as from serendipitous obser- 
vations close to the Galactic Centre, obtained in 1996 (net 
observing time of 851ksec), 1997 (785ksec) and north- 
ern spring 1998 (414ksec). Data covering (part of) the 
Galactic Centre are routinely checked for X-ray bursts. 
We found a total of 63 bursts from KS 1731—260 in the 
above interval. 

The Rossi X-ray Timing Explorer (RXTE) was 
launched 1995 December 30. This satellite carries three 
instruments, including the All-Sky Monitor (ASM), which 
consists of three Scanning Shadow Cameras (SSCs) 
mounted on a rotating drive such that one is perpendic- 
ular to the other two (Levine et al. 1996). The cameras 
are held stationary for 90-s intervals, called 'dwells', dur- 
ing which data are accumulated. The mount rotates the 
assembly through a 6° angle between dwells until a full 
rewind is necessary. In this manner, ~80% of the sky is 
observed every 90-min orbit. Each dwell provides intensi- 
ties in the 1.5-12 keV band for all known sources in the 
field of view of each camera, and results that meet a set 
of reliability criteria are saved in electronic tables that are 
available over the world wide web[]. An alternate, time- 
series, data mode records the total counts in each SSC in 
0.125-s bins, but does not preserve the position of each 
photon (Levine et al. 1996, Smith et al. 2001). Over the 
~5.5 years of RXTE operation, KS 1731-260 was observed 
an average of ~10 times a day, yielding ~17000 individual 
intensity measurements. 

4. Analysis and results 

4.1. A long X-ray flare 

Following the discovery of the long X-ray flare in 
4U 1735-44 (Cornelisse et al. 2000), we found two other 
such events in our WFC database of X-ray bursters: one 
from Ser X-l, presented by Cornelisse et al. (2001), the 
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Fig. 1. * 

a. RXTE/ASM (1.5-12 keV) light curve of KS 1731-260 from 
1996 January 6 (MJD 50088) to 2001 May 1 (MJD 52030). 
The data shown are two-day averages and normalized to 
the count rate of the Crab (75ctss -1 SSC -1 ); only intervals 
with more than one dwell per day were included, b. The 
BeppoSAX/WFC (2-28 keV) data during the interval 1996 
August 12 to November 1 (MJD 50307-50388; marked by the 
dotted lines in panel a). Data points (typically averages per 
satellite orbit) which were less than half a day apart are 
connected to guide the eye, and are normalized to the Crab 
(2ctss _1 cm -2 ). The gray data points are the one-day aver- 
ages (with more than one dwell per day) of the RXTE/ASM 
in the same interval. A long X-ray flare in the WFC data is 
seen just before MJD 50350 (1996 September 24). The vertical 
tickmarks indicate the times of the 20 X-ray bursts seen with 
the WFC. They all occur before the flare, c-e. WFC light 
curves around the flare (marked by dotted lines in panel b) for 
channels 1-11 (2-5 keV; c) and 12-31 (5-28 keV; d), and the 
hardness ratio of the two above defined energy bands (e) . The 
time resolution is 5min. Overplotted in c with open squares 
are the data from individual RXTE/ASM dwells. 



other from KS 1731—260, presented here. The long X- 
ray flare of KS 1731-260 can be seen clearly in Fig. [l|b, 
which shows the WFC data obtained on KS 1731—260 in 
1996, together with the one-day averages of the ASM. A 
blow-up of this light curve in two energy bands near the 
flare is shown in Figs. Qc and d. The flare started near 
MJD 50349.42 (1996 September 23.42 UT), the observed 
maximum was seen about half an hour later. Subsequently 
the source decayed back exponentially to its pre-flare level 
with a decay time of r cxp ~ 2.7hr (2-28 keV; see Table 1). 
The decay is faster at high energies than at low energies 
(Table 1). The distinction in the two decay light curves 
is shown clearly in the hardness curve (Fig. ||e): the spec- 
tra harden during the rise to maximum, and soften again 
during the subsequent decay. 

The flare was partly covered by the RXTE/ASM 
(Fig. 0c) and corresponds to the highest data point in 
the ASM light curve (Fig. |l]a). This enhancement was al- 
ready noticed by Wijnands (2001), who, however, could 
not determine its nature because no evidence for soften- 
ing could be found. We note that there are other high data 
points which stand out in the ASM light curve (e.g. near 
MJD 50623: 1997 June 24), which might be remnants of 
other long X-ray flares. We examined the time-series data 
mode for each dwell that yielded an anomalously high in- 
tensity for KS 1731—260 and found no evidence that a X- 
ray burst had occurred during any of these observations, 
which would increase the mean count rate during a dwell. 
Although some of the high excursions occur when the Sun 
was close in position in the sky to KS 1731—260 (e.g., 
near MJD 50400: 1996 November 13), suggesting that the 
increase is not due to the source itself, we can not rule out 
that the other high data points were obtained during long 
X-ray flares (if so, this would indicate recurrence times of 
less than a year). 

In Fig. ^a, we show the count rates (2-28 keV) of the 
observation near the start of the X-ray flare. The data 
show a short increase in intensity for about 4 s, about 380 s 
after the start of the observation. We will refer to this as 
'a precursor' hereafter. Afterwards, the intensity decreases 
somewhat erratically to almost its pre-burst level^J. Just 
before the data gap, the flux may start to rise again, but 
unfortunately we cannot be certain. In Fig. [|a, we again 
show the count rates, now including the X-ray flare, on a 
logarithmic time scale. After the gap (~13min later), the 
flux is high and it is possible the flare has already passed 
its maximum. 

We fit the decay light curve of the precursor with an ex- 
ponential and derive the peak intensity, integrated count 
rate, and exponential decay time, r cxp . In Table 1, we 
compare these numbers with what is found for the other 
63 X-ray bursts seen with the WFC at comparable persis- 
tent pre-burst fluxes (i.e. before MJD 51000). We find that 
the characteristics of the precursor are comparable to the 



2 In Figs, [yc and d, the flux appears to increase steadily 
during the start of the flare. This is due to the choice of data 
rebinning. 
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Fig. 2. a. WFC (2-28 keV) count rates of the observation 
just before the long X-ray flare. Time = s corresponds 
to UT 1996 September 23 09:52:02 (MJD 50349.41113); 
the time resolution is 2 s. Clearly, a short and weak X- 
ray burst, which we call a precursor, can be seen 380 s 
after the start of the observation. The arrow points to 
the zero time scale of Fig. 3. b. Example of a WFC (2- 
28keV) observation around another short and weak X- 
ray burst. Time = 0s corresponds to UT 1997 April 13 
19:15:00 (MJD 50551.80208); the time resolution is 2s. 



other bursts, but note that the data are not of sufficient 
quality to determine whether there occurs significant soft- 
ening during its decay (thus, we cannot prove it is a type I 
burst). For comparison, we show a light curve around an- 
other similar short and weak X-ray burst (one of the 63 
bursts) in Fig. ||b, to illustrate that such bursts are not 
uncommon in KS 1731—260 (see also Muno et al. 2000). 

In Fig. |l]b, we show with tickmarks the times of oc- 
currence of normal X-ray bursts observed with the WFC 
around the time of the long X-ray flare (no bursts were 
found during the short ASM dwells on the KS 1731-260 
region in the same interval). We found 20 normal bursts, 
all of which occur before the flare. The mean burst rate, 
corrected for data gaps, is ~3 per day, but the wait- 
ing times between bursts are variable, the shortest being 
-0.08 day. 

No normal bursts were found during the WFC obser- 
vations which were obtained up to ~35 days after the flare, 
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Fig. 3. a. WFC (2-28 keV) count rates of the large flare 
using a logarithmic time scale, with the data rebinned log- 
arithmically in time. Time = s corresponds to UT 1996 
September 23 09:58:21 (MJD 50349.41552). The dotted 
line shows the mean count rate before the large X-ray flare, 
b— e. Results of black-body fits to the net burst emission 
(= total burst emission minus pre-burst emission): bolo- 
metric black-body flux, Ft,b, in units of 10 -8 ergcm" 2 s" 1 
(b), black-body temperature, fcTbb (c), apparent black- 
body radius at 7kpc, Rhb,i (d), and the goodness of fit 
expressed in reduced \ 2 f° r 23 degrees of freedom (e) . 



despite considerable exposure time (227ksec). Given that 
20 bursts occurred during 589 ksec of observation time be- 
fore the flare we would have expected to find ~8 bursts 
after the flare. We verified the significance of finding no 
bursts after the flare by performing Monte Carlo simula- 
tions. We randomly varied the burst waiting times be- 
tween the lowest observed waiting time (0.08 day) and 
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0.6 day (i.e. symmetric around 0.34 day). We then deter- 
mined the number of bursts which occurred within the 
observing windows after the flare. By doing 10 6 simula- 
tions wc found that the expected number of bursts is in- 
deed 7.6±2.4 (1 ex); the probability of observing no bursts 
after the flare is lxlO -4 . We verified our method by ap- 
plying the same procedure for the interval before the flare: 
20±3.5 are expected; indeed 20 are seen. 

Another way to estimate the significance is to as- 
sume that the number of observed normal bursts follows 
a Poisson distribution. In this case we find that the prob- 
ability of seeing no bursts instead of the 8 expected is 
3.4xl0~ 4 . Finally in WFC observations during 1997 (first 
observation was on MJD 50509, i.e. 1997 March 2) and 
1998, when the source was at comparable flux levels, the 
source was also bursting, again with similar burst frequen- 
cies as those before the flare (see Fig. ||). 

Given the above, we conclude that the burst activity 
was significantly reduced or had ceased after the long X- 
ray flare. 



We estimated the fluence, E^, of the flare by assuming 
that the rise from the start to the observed maximum is 
linear and that the decay after this maximum is described 
by an exponential. In Table 1, we list the values of E^, the 
maximum observed net-burst black-body flux, -Fbb,max ? 
the ratio r = E h /F hhtIaax , and the ratio 7 = F pcrs /F b b,max 
(where F pers is the bolometric persistent source flux ap- 
proximated by extrapolating the observed persistent spec- 
trum between 0.01 and 100 keV). Note that the values of 
r and 7 are rather uncertain since we missed the peak 
of the flare. We find that the fluence E\> of the flare is a 
factor of 250-425 larger than what is observed for normal 
X-ray bursts in this source (see e.g. Muno et al. 2000). 
We note that if we do not use the exponential but simply 
integrate the observed fluxes during the flare and inter- 
polate over gaps, we obtain a slightly higher fluence of 
~ 2.4 x 10~ 4 erg. The flare does not reach the Eddington 
limit during our observations: fbb.max — 0.4 L^dd given 
the Eddington flux of 6.3 x 10~ 8 erg cm' 2 s' 1 inferred 
by Muno et al. (2000) from radius expansion bursts of 
KS 1731-260. 



4.2. Spectral behaviour 

In order to derive the characteristics of the long X-ray flare 
emission, we created spectra throughout its duration, in- 
cluding the precursor. The time intervals were chosen such 
that we had sufficient statistics per spectrum. The net flare 
spectra (i.e. total source spectra minus pre-flare persistent 
emission spectrum C described below) were satisfactorily 
fitted with emission from a black-body (x 2 ed ^ f° r ^3 
degrees of freedom) . We note that some of the x? c d va l ues 
are rather small, perhaps due to a slight overestimate of 
the systematic errors. As we cannot constrain interstellar 
absorption well with the WFC, we fixed Nn to the value 
derived by Barret et al. (2000), 1.3 xlO 22 atoms cm" 2 . The 
results are shown in Figs. ^|b-e. The black-body temper- 
ature reaches a maximum of kT — 2.4 ± 0.1 keV and de- 
creases thereafter. The apparant black-body radius at the 
distance of 7 kpc remains more or less constant through- 
out the burst with a value of ~6 km; as expected, the re- 
lation between kT and Fbb is consistent with a power law, 
FbbOcT", with a = 3.9±0.4. Similar black-body parame- 
ters are found for normal X-ray bursts from this source 
(e.g. Muno et al. 2000). 

We also performed spectral fits to the flare spectra 
not subtracting the pre-flare persistent emission. The fits 
to the spectra up to 12200 s into the flare were accept- 
able: Xrcd = °- 7 - L2 (maximum F bb = 2.9 ± 0.1 x 
10~ 8 ergcm~ 2 s _1 , maximum kT = 2.2±0.1keV, i? bb ,7 ^ 
10 km). However, the fits to the spectra of the remainder 
of the burst were not acceptable: x 2 cd = 2.8 — 7.7. This 
indicates that while it seems more likely the persistent 
emission was present during the flare, we can not rule out 
its absence during the first ~3.5hr. In the remainder of 
the paper we use the result from the spectral fits to the 
net-flare emission. 



4.3. Persistent emission 

The change in bursting behaviour just after the flare led 
us to investigate whether the persistent emission is influ- 
enced as well. We calculated hardness ratios (as defined in 
Sect. 3.1) during the interval 1996-1998, when the source 
was at comparable persistent flux levels. These are shown 
in Fig. ^; clearly, the source did not have exceptional hard- 
ness values after the flare. 

In order to quantify the above, we created a total of six 
X-ray spectra of the persistent emission before and after 
the flare during different intervals between MJD 50308- 
50388 (see Table 2). For each spectrum data from vari- 
ous contiguous observations with the same pointing were 
accumulated. To parametrize the persistent emission we 
fitted the 2-20 keV spectra to an absorbed cut-off power- 
law model plus a gaussian line at 6.4 keV (see Barret 
et al. 2000). As above, we fixed the value of Nh to 
1.3 xlO 22 atoms cm' 2 ; furthermore, we fixed the line width 
of the gaussian, a, to 0.1. This provided a good descrip- 
tion of the spectra^; the results are presented in Table 2. 
The unabsorbed persistent luminosity before and after the 
flare was ~ 0.1L E dd (2-10 keV). 

We conclude that the persistent X-ray spectra days to 
weeks after the flare (when no normal bursts occurred) 
are not different from the spectra days to weeks before 
and months to years after the flare (when normal bursts 
did occur). 

5. Discussion 



3 Previously, a black-body component (or a disk black-body 
component) has been included in spectral fits (Barret et al. 
2000; Narita et al. 2001). We, however, found their inclusion 
made the fits rather unstable, while providing little (spectra A 
and F) or no (B-E) improvement. 
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Table 1. KS 1731—260 precursor and superburst characteristics 





Precursor 


Type I bursts 


peak rate" 


2.3+0.4 ctss" 1 cm" 2 


1.8-5.3 (3.6) ctss^cm" 2 


integrated 






count rate 


~13cts cm" 2 


13-45 (29) ctscnr 2 


r oxp (2-28 keV) 


5.6± 2 : 8 s 


2-11 (5.5) s 



Superburst 



parameter 


value 


Xrod/dof 


r cxp (2-28 keV) 


2.7±0.1hr 


0.94/116 


(2-5 keV) 


4.4±0.3hr 


0.7/116 


(5-28 keV) 


2.1+0.1 hr 


0.8/116 


parameter 


at earth 


at 7kpc 


F b 


1.0+0.1 x 10~ 8 erg s _1 cm -2 


5.9±0.8xl0 37 ergs" 1 


E b 


~1.7xl0~ 4 erg cm~ 2 


~9.8xl0 41 erg 


Fbb,max 


2.4±0.1 x 10~ 8 erg s _1 cm" 2 


1.4±0.1xl0 38 ergs" 1 


parameter 


value 




T=Eb/Fbb,max 


~2.0hr 




'Y — Fpcrs/Fbb,max 


0.41±0.06 





a For the normal bursts we give the range, and the average value in parentheses. 
b Unabsorbed pre-superburst luminosity (0.01-100 keV). 



Table 2. Persistent emission spectral fit (absorbed cut-off power- law) parameters" 





A 


B 


C 


D 


E 


F 


Interval (MJD-50000) 


318.6-322.2 


323.4-326.0 


348.7-349.3 


350.2-250.5 


366.7-369.0 


383.1-385.0 


i^-iokev (10~ 9 ergs -1 cm" 2 ) 


4.6±0.2 


4.8±0.3 


5.2±1.1 


5.7±1.4 


4.3±0.4 


4.2±0.8 


r b 


1.04±0.03 


1.20±0.05 


1.0+0.1 


1.2±0.2 


1.03±0.06 


0.9±0.1 


E cut (keV) c 


4.3±0.1 


4.6±0.2 


5.2±0.7 


5.6±1.0 


4.5±0.2 


3.4±0.3 


norm. d 


1.34±0.03 


1.63±0.05 


1.3±0.1 


1.5±0.2 


1.17±0.05 


1.26±0.09 


gnorm. e (xl0~ 3 ) 


4.1±0.8 


6.6±1.3 


_/ 


_/ 


_/ 


_/ 


X 2 /dof 


1.6/21 


0.9/21 


1.2/22 


1.1/22 


1.0/22 


1.5/22 



° Nh was fixed to 1.3x10 atoms cm , see text. 
6 Power-law photon index. 
c Cut-off energy. 

d Power-law nomalisation in photons cm~ 2 s _1 at 1 keV. 

' Normalisation in photons cm -2 s _1 in the line modeled by a Gaussian fixed at 6.4 keV and width (a) of 0.1 keV. 



* Including a gaussian did not improve the fit significantly. 



5.1. Very long thermonuclear runaway events 

In many of its characteristics, the long X-ray flare resem- 
bles type I bursts from KS 1731—260 (e.g. Muno et al. 
2000) as well as type I bursts from other sources. The 
rise to maximum is relatively fast with respect to the 
decay, the decay is exponential, and the spectra soften 
(indicating cooling) during the decay. This strongly sug- 
gests a thermo-nuclear runaway event on the neutron star. 
The duration of the X-ray flare, however, is about half a 
day, i.e., a factor of ^1500-4000 longer than normal X-ray 
bursts from this source. Therefore, the fluence, ~ 10 42 erg 
(at 7kpc), is also much larger than typical. The longest 
'normal' X-ray bursts seen in general are much shorter 
and less energetic, and all are radius expansion events: 
they had durations on the order of half an hour (Hoffman 
et al. 1987: probably from 4U 1708-23; Kuulkers et al. 
2001: GX 17+2), fluenccs of ~10 41 erg, exponential decay 
times of 5.5 min (4U 1708-23) and 4.6 min (GX 17+2), 
and ratios of the fluence to peak net-burst flux (r) of 



~5.4min and 6.6 min, respectively. These values are an 
order of magnitude smaller than those for the superbursts. 

The long X-ray flare resembles the superbursts seen 
in five other sources (4U 1820-30, Scr X-l, 4U 1735-44, 
4U 1636—53 and GX 3+1), as can be seen in Table 3, 
where we compare the properties of all seven flares ob- 
served to date. The duration of the superbursts are es- 
timated by eye from the light curves and the hardness 
curves. With respect to duration of about half a day the 
flare in KS 1731—260 is the longest among the superbursts 
seen so far. 

The superburst seen in 4U 1820-30 (Strohmayer 2000; 
Strohmayer & Brown 2001) is the only one for which the 
full rise to maximum has been observed. That superburst 
was immediately preceeded by a normal X-ray burst. This 
'precursor' lasted for ~20s. Both the 'precursor' and the 
superburst exhibited radius expansion, indicating the lu- 
minosity reached the Eddington limit. The superbursts 
from all sources but 4U 1820—30 seem to be a factor 2- 
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Table 3. Properties of superbursts ordered along their exponential decay time a 



source 


4U 1820 30 


Ser X-l 


4U 1735 44 


4U 1636 53 


GX 3+1 


KS 1731 260 


instrument 


PCA 


WFC 


WFC 


ASM,PCA 


ASM 


WFC,ASM 


energy range 


2-60 keV 


2-28 keV 


2-28 keV 


1.5-12 keV 


1.5-12 keV 


2-28 keV 


duration (hr) 


>2.5 


~4 


~7 


£1-3 


>3.3 


~12 


precursor? 


yes 


? 


? 


? 


? 


yes 


"T eX p (hr) 


~1 


1.2±0.1 


1.4±0.1 


1.5±0.1/3.1±0.5 


1.6±0.2 


2.7±0.1 


ipors (LEdd) 6 


~0.1 


~0.2 


~0.25 


-0.1 


~0.2 


~0.1 


fcT max (keV) 


~3.0 


2.6±0.2 


~2.6±0.2 


? 


~2 


2.4±0.1 


Lpcak (10 38 ergs) c ' d 


3.4±0.1 


1.6±0.2 


1.5±0.1 


~1.2 


~0.8 


1.4±0.1 


E h (10 42 erg) 


£1.4 


~0.8 


£0.5 


~0.5-l 


£0.6 


~1.0 


r=E b /L peak (hr) d 


£1.1 


~1.4 


£0.9 


-1.2-2.3 


£2.1 


~2.0 


T— Lpers /Lp ea k 


~0.1 


~0.4 


~0.4 


~0.2 


~0.5 


~0.4 


tno bursts (days) e 


7 


~34 


>7.5 


? 


? 


>35 


H/He or He donor 


He 


? 


H/He 


H/He 


? 


? 


references-^ 


S00,SB01 


C01 


COO 


W01 


K01 


this paper 



a A question mark denotes an unknown value. 

b We used the 0.01-100 keV unabsorbed flux from spectral fits and the observed maximum flux 

during radius expansion bursts. 
c Unabsorbed bolometric peak (black-body) luminosity. 

d Since the rise to maximum was only observed for 4U 1820—30, the values for the other sources are to be used with caution. 
e Time of cessation of normal X-ray bursts after the superburst. 

/ C00=Cornelisse et al. 2000, C01=Cornelisse et al. 2001, K01=Kuulkers 2001, S00=Strohmayer 2000, 
S01=Strohmayer & Brown 2001, W01=Wijnands 2001. 



CD 
£Z 
XJ 

O 



CD 
O 



•i 



400 600 800 

Days (MJD - 50000) 

Fig. 4. WFC hardness (as defined in Fig. le) curve from 
MJD 50275-50950 (1996 July 11 - 1998 May 17) when 
the persistent X-ray intensity was at comparable levels. 
Data points are 1-day averages. The dotted line denotes 
the time of the long X-ray flare. In the upper part we 
have marked the times of the bursts with bars. Clearly, no 
bursts are found just after the flare. Normal bursting had 
resumed in the next observation. 



4 fainter than their brightest normal bursts, although we 
should stress again that only in 4U 1820—30 the rise to 
maximum was fully covered. Next best is KS 1731—260, 
for which we observed the start and found a short (~4 s) 
and weak burst, most probably a type I burst, preceding 
the superburst by ^200 s. Note that the precursor does 
not necessarily need to be a trigger of the superburst (see 
also Strohmayer & Brown 2001). The observation of pre- 



cursors close to the superbursts is another difference from 
normal (type I) bursts, which do not have such precursors 
(see e.g. Lewin et al. 1993). 

We found evidence that the occurrence of normal 
bursts stopped for more than a month after the super- 
burst in KS 1731—260. Similarly, no bursts were observed 
for about a week after the superburst in 4U 1735—44, 
whereas in Ser X-l no bursts were seen for ^5 weeks 
after the superburst (Cornelisse et al. 2000, 2001; see 
also Table 3). This suggests that the superburst has pro- 
found influence on the outer envelope of the neutron star 
where normal bursts originate. The outer envelope may 
be heated by the energy generated by superburst, causing 
the burning of hydrogen and/or helium to become stable 
(see Cumming & Bildsten 2001). However, the persistent 
emission of KS 1731—260 during the non-bursting interval 
is not significantly different compared to intervals of nor- 
mal bursting. This seems evidence against the idea that 
the neutron star is significantly hotter, but might also just 
be because the observed persistent emission is dominated 
by emission from the accretion disk. 

5.2. Unstable carbon burning... 

Assuming unstable nuclear processes, we can use the time 
scales and total energy to estimate where the burning oc- 
curs, how much matter needs to be involved, etc. (see also 
Cumming & Bildsten 2001; Strohmayer & Brown 2001). 
Since the fusion of elements higher than hydrogen releases 
about lMeV per nucleon (~ 10 18 ergg _1 ; see, e.g., Lamb 
& Lamb 1978), one needs to process roughly 10 24 g of 
matter to produce the total energy of 10 42 erg released 
during a superburst. With the mass accretion rates of 
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~10" 9 M© yr _1 (as appropriate for sources which radiate 
at ~0.1-0.3 iEdd) the expected recurrence time is at least 
half a year (depending on how much of the accreted mat- 
ter is consumed during a superburst). This indicates that 
they are indeed rare events. 

The above mass estimate corresponds to an ignition 
column depth y = AM/47ri? 2 s ~ 8 x 10 10 gcm~ 2 (where 
AM is the amount of burned matter and R ns the radius 
of the neutron star; we take R ns — 10 km). At these col- 
umn depths, the thermal time is on the order of 1.5-3 hrs, 
almost independent of the temperature (see e.g. Brown 
& Bildsten 1998). This is in excellent agreement with the 
decay time scale of the superbursts. We note that the den- 
sity at column depths ~ 8 x 10 10 g cm~ 2 is p ~ 10 8 g cm~ 3 
(see, e.g., Brown & Bildsten 1998). 

The above arguments have led to the suggestion that 
superbursts are due to unstable carbon burning deeper in 
the neutron star, in a layer beneath the (un)stable hydro- 
gen/helium or helium layer (Cumming & Bildsten 2001; 
Strohmayer & Brown 2001). Obviously, in order for this 
to work, the ashes of the burning hydrogen/helium layer 
need to contain carbon. One way to achieve this is to have 
stable burning of helium. This may occur in 4U 1820—30, 
in which the neutron star is thought to accrete pure he- 
lium because the companion is most likely a helium white 
dwarf (based on the very short orbital period (Rappaport 
et al. 1987), which is consistent with the short bursts this 
source shows. 4U 1820—30 shows long intervals in which 
no normal X-ray bursts occur, which gave ample time to 
generate the necessary amount of carbon (Strohmayer & 
Brown 2001). 

The above mechanism cannot work for all systems, 
however. In particular, for 4U 1735-44 and 4U 1636-53 
it has been shown that the donors provide hydrogen-rich 
material (Augusteijn et al. 1998). For burning of a hydro- 
gen/helium mix, one expects only small amounts of carbon 
to be produced, both for the case that the burning is sta- 
ble and for the case that it is unstable (e.g. Schatz et al. 
1999, 2001). In particular, during normal bursts temper- 
atures in excess of ~ 10 9 K are reached and a break-out 
of the CNO cycle occurs, creating elements far beyond 
the iron group. Nevertheless, over a sufficiently long time 
(many years for mass accretion rates of a few tenth of 
Eddington) sufficient carbon may be produced to cause a 
superburst (Cumming & Bildsten 2001). 

Another issue is whether the burning of pure car- 
bon can be unstable at such low column depths as 
8xl0 10 gcm~ 2 . Brown & Bildsten (1998) found ignition 
at >5xl0 12 gcm~ 2 for pure carbon at Eddington accre- 
tion rates. The superbursters accrete near a tenth of the 
Eddington accretion rate, making it even more of a chal- 
lenge to ignite carbon at the lower pressures implied by the 
observations. Strohmayer & Brown (2001) show that one 
may underestimate the amount of energy released when 
integrating the superburst X-ray flux, and therefore also 
the amount of carbon consumed during the superburst. 
According to the authors, most of the energy is released as 
neutrinos, and most of the heat produced by the burning 



largely flows inwards. If so, the large amounts of carbon 
required (i;10 26 g), imply very long recurrence times for 
the superbursts in 4U 1820— 40 (decades). An alternative 
suggestion is that the precursor helium flash triggered un- 
stable carbon burning, but Strohmayer & Brown (2001) 
consider this to be less likely. 

Cumming & Bildsten (2001) have a different view. 
They point out that the heavy elements in the ashes will 
reduce the thermal conductivity, which may lead to large 
temperature gradients and ignition at the observed col- 
umn depth. They find this works for mass accretion rates 
above 10% of the Eddington rate and for mass fractions of 
carbon larger than 0.1. They argue that the fact that up 
to now superbursts have only been reported in systems 
with accretion rates between 10 and 30% of Eddington, 
is that at lower mass accretion rates the superbursts are 
easier to distinguish from the persistent emission. 

Cumming & Bildsten (2001; see also Strohmayer & 
Brown 2001) argued that the cooling following the car- 
bon burning may provide a large heat flux into the hydro- 
gen/helium burning layers, stabilizing the burning there. 
For mass accretion rates of 10-20% of Eddington bursting 
would then cease for about 15-20 days after the super- 
burst. This is more or less consistent with what is observed 
in KS 1731-260 and ScrX-1. 

5.3. ... or runaway electron capture? 

Although with unstable carbon burning it appears possi- 
ble to understand the superbursts, there may be an in- 
teresting alternative. The densities at which the Fermi 
energy of the electrons, which provide the pressure in 
the outer layers, becomes equal to the difference in rest- 
mass of a neutron and a proton are 2-3xl0 7 gcm~ 3 
(where the exact number depends on the composition). 
This density corresponds to a column density of y ~ 
2 x 10 10 57 4 1 3 gcm~ 2 , where 314.3 is the surface gravity nor- 
malised to 10 14 ' 3 cms -2 , appropriate for a neutron star 
mass M ns = 1.4 M Q and R ns = 10 km (e.g. Bildsten & 
Cumming 1998). This is not too different from the values 
inferred above. Hence, if any protons are present, they will 
capture electrons via the p(e~ ,v)n reaction (Rosenbluth 
et al. 1973; see also Taam et al. 1996; Bildsten & Cumming 
1998). The released neutrons are then captured by heavy 
nuclei (rather than by protons, see also below) releas- 
ing ~7-8MeV per nucleon (Bildsten & Cumming 1998). 
Depending on the temperature (which sets the precise dis- 
tribution of electron energies) and the thermal conductiv- 
ity, this process may be unstable: an increase in temper- 
ature leads to a larger fraction of the electrons at high 
energies; if the region does not cool sufficiently fast, this 
leads to more captures, a further increase in temperature, 
etcetera. 

There have been some studies of the above described 
mechanism, but only for the case in which no heavy nu- 
clei are present and the neutrons are captured by protons. 
This process is called 'deep hydrogen burning' (Hameury 
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et al. 1982, 1983; Woosley & Weaver 1984; Fushiki et al. 
1992; Taam et al. 1993, 1996). This type of burning can 
become unstable when the neutron star envelope has a 
low temperature (^10 7 K) and/or when the abundance of 
CNO elements is very low (Z <, 0.001). The latter may 
not be as unrealistic as it appears at first, since the origi- 
nal CNO abundance received from the donor may be sig- 
nificantly reduced due to nuclear spallation processes in 
the accretion shock (Bildsten et al. 1992). In the mod- 
els, long (>30min) X-ray bursts have been produced at 
moderate accretion rates (~0.1 MEdd), with total nuclear 
energy releases of up to several times 10 41 erg (Woosley 
& Weaver 1984, their model 6; Fushiki et al. 1992). This 
is still an order of magnitude less than observed in super- 
bursts (but rather similar to the 'long' X-ray bursts seen 
from 4U 1708-23 and GX 17+2!). 

Several other questions have to be answered as well, 
however. Foremost, whether there is (enough) hydrogen 
left after the (un)stable hydrogen/helium burning. Various 
calculations suggest residual hydrogen mass fractions of 
Xh ~0.04-0.2 for mass accretion rates of >1O _9 M0 yr _1 
(e.g. Taam et al. 1996; see also Koike et al. 1999). (Of 
course, this mechanism cannot explain the superburst in 
4U 1820—30, since for that source there are no protons to 
begin with. We note, however, that this superburst seems 
to be somewhat different from the others, so a different 
mechanism may not be unexpected.) 

To produce 10 42 erg one needs ~10 23 g of residual hy- 
drogen (given the energy release of 7-8 MeV per nucleon). 
Thus a total accreted matter of 5-25 x 10 23 g is needed; for 
accretion rates of ~10~ 9 M Q yr~ 1 the recurrence times of 
superbursts would then be 0.25-1.25 yr. This is in an or- 
der of magnitude shorter than the case for unstable carbon 
burning as estimated by Cumming & Bildsten (2001). We 
estimate the maximum expected energy release for run- 
away electron capture from E max = ATrR^ s yE nuc XH = 
2 x 10 41 i? 4 s 10 Xh,o.i erg, where R ns — 10i? nSj iok m is the 
radius of the neutron star, E nuc the above quoted energy 
release, and Xh = O.lVyo.i (the fourth-power depen- 
dence on radius arises because of the dependence of y on 
the surface gravity). This may be consistent with the ob- 
served value of 10 42 erg, if the radius of the neutron star 
is larger than 10 km, as is indeed inferred for many equa- 
tions of state (see, e.g., the review by Lattimer & Prakash 
2000), and/or the residual mass fraction of hydrogen is 
larger than 0.1. Thus, unstable electron capture appears 
a viable alternative mechanism for the superbursts. 
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